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Abstract 

A microscopic formulation of the definition of both the heat flux and the viscous stress 
tensor is proposed in the framework of kinetic theory for relativistic gases emphasizing on 
the physical nature of such fluxes. A Lorentz transformation is introduced as the link be- 
tween the laboratory and local comoving frames and thus between molecular and chaotic 
velocities. With such transformation, the dissipative effects can be identified as the aver- 
ages of the chaotic kinetic energy and the momentum flux out of equilibrium, respectively. 
Within this framework, a kinetic foundation of the ensuing transport equations for the 
relativistic gas is achieved. To our knowledge, this result is completely novel. 



I. INTRODUCTION 

Relativistic kinetic theory is not a new subject, however it now finds itself in 
a spothght due to the increasing interest in relativistic thermodynamics triggered 
by recent heavy ion collisions experiments, electron-positron plasma generation, and 
the traditional astrophysical applications of relativistic hydrodynamics. The theory 
has its roots in original works by Jiittner [ij for the equilibrium case while the first 
kinetic theory treatment was formulated by Israel |2|. In such work, the author 
finds an expression for the stress energy tensor by solving the Boltzmann equation 
using a Chapman- Enskog expansion. However, since the systematic (hydrodynamic) 
and chaotic (or peculiar) components of the total velocity of a given molecule are not 
explicitly distinguished, the different contributions to this tensor cannot be identified 
in the same fashion as in the non-relativistic case (3|, \^. Instead, projections in 
parallel and orthogonal directions with respect to the hydrodynamic velocity of the 
stress energy tensor are used and the interpretations of the different contributions 
agree with those that follow from the phenomenological counterpart as developed by 
Eckart |5|. This procedure is essentially followed by most authors 6|, l7|. 

On the other hand, in non-relativistic kinetic theory, a clear distinction can be 
made between the effects caused by the "bulk", or mechanical, properties of the fluid 
and its microscopic ones. This permits the identiflcation of dissipative fluxes, i. 
e. heat and viscosity effects, as averages of chaotic quantities jj|. In particular, the 



interpretation of heat flux as the average o 



the chaotic kinetic energy flux, as deflned 



more tha,> a century ago by R^ C.ausius ^ -d J, C, Maxwell H. .s assorted, Th,s 

concept is absent in the relativistic case, as was clearly noted in Ref. |11|. In that 
work, a flrst proposal of how to introduce the chaotic velocity concept in relativistic 
kinetic theory was put forward. In this work, we follow the same line of thought and 
take it a step forward by explicitly introducing Lorentz transformations in the stress- 



energy tensor integral in order to separate mechanical and chaotic effects. By doing 
so, we are able to clearly define the heat flux and the viscous stress tensor as the 
average of chaotic energy and momentum fluxes in an arbitrary frame, respectively. 
To accomplish this task we have divided this work as follows. In Section II, we 
briefly review the non-relativistic setup for calculating the dissipative fluxes. The 
relativistic framework is introduced in Section III where Lorentz transformations 
are used in order to introduce the chaotic velocity and obtain the corresponding 
expressions for the heat flux and viscous stress tensor. A brief discussion of the 
results and flnal remarks are included in Section IV. 

II. NON-RELATIVISTIC KINETIC THEORY 

Kinetic theory serves as the microscopic foundation of irreversible thermodynam- 
ics and is capable of producing both the system of transport equations as well as 
the constitutive equations needed in order to make it a complete set describing the 
dynamics of fluids J3|, 1^. As usual, the distribution function / {f,v,t) is such that 
/ (r, V, t) dfdv is the number of molecules contained in a 6-box in the phase space 
corresponding to position r and molecular velocity v. The local variables are thus de- 
flned as averages weighted by this function. The local particle density, hydrodynamic 
velocity and energy density are thus deflned as 

n = I f^^^d^v (1) 

^ ' 7/(o)d\ (2) 



u = — I V 
n 

e = - [ -mv'^f^'^d'v (3) 

n I 2 



respectively, where f^^^ is the local equilibrium distribution function: 

T being the temperature, m the molecular mass and ks the Boltzmann constant. 
The evolution of the distribution function is given by the Boltzmann equation. For 
a simple (one component), non-degenerate, diluted gas in the absence of external 
fields the kinetic equation reads 

| = i(//') (5) 

where, if g and a are the relative velocity and cross section for a collision between 
two particles respectively, the collision term is given by 

J (/, f) = j j if '/i ' - / h}9<ydndvl (6) 

Primes denote quantities after the interaction and Vt is the solid angle. The well 
known Maxwell-Boltzmann distribution function given in Eq. (j4]) is precisely the 
solution of J (/, /') = 0, namely the homogeneous Boltzmann equation. The solu- 
tion to the inhomogeneous, out of equilibrium, case is in general obtained via the 
Chapman-Enskog method in which the general solution is written as 

/ = /(°) + f'^ (7) 

where the second term contains corrections to the equilibrium solution to first order 
in the gradients of the local variables and gives rise to the dissipative fluxes. This 
term includes only dissipative effects once the solubility constraints 

/Wrf\= [vfi^U'v= /t;2/«rf\ = (8) 



are introduced such that the local variables are defined through the local equilibrium 
distribution solely. 



In this framework, the transport equations are obtained by multiplying Eq. ([5]) 
by a collision invariant and integrating over v. Such procedure yields the Maxwell- 
Enskog transport equation 

— f il)f£v + V ■ /" i)vf£v = (9) 

which accounts for particle, momentum and energy balances for ip = 1, v, v"^, respec- 
tively. Indeed, taking ip = 1 in Eq. ([7]) yields the continuity equation 

_ + VM = (10) 

For ip = V one obtains 

^ + V.r = (11) 

ot 

where we introduced the stress tensor 



^^/ 



vvfdh (12) 



Finally, the energy balance is obtained for if) = v^: 

One 



+ V ■ Je = (13) 



^'h. 



dt 
where we have defined the total energy flux as Jp = f v'^vfd^ 

poses the molecular velocity in its two basic components, usually written as 

v = u + k (14) 

where k is the chaotic or peculiar component. In this case such expression arises in a 
very natural way by observing the argument of the exponential function in Eq. (^. 
It is clear that J kfd^v = in view of Eqs. (jl]) and (ITO|) and thus 

e=-u^ + e (15) 



7 = nuu + nksTl 



+ v 



Je = -nv^u + nue + nu ■ {nksTl + t ) + q 



(16) 
(17) 



where 



ne 



is the internal energy density per particle and the dissipative fluxes are given by 



^-> 



q 



kkf^^^d^k 



-kf^^d^k 



(19) 



(20) 



Introducing these definitions in the transport equations and using the local equilib- 
rium assumption, one obtains the well known set of hydrodynamic equations for the 
non-relativistic fiuid. 

In the next section, it will be shown how these ideas can be extrapolated in a 
very natural way to the relativistic framework. In order to make the transition more 
clear we want to point out at this stage the key role of the transformation v = u + k 
in the formalism. Notice that such a transformation can also be expressed in terms 
of a Galilean matrix in space-time, that is 



where the Galilean transformation is given by 

/ 1 u^/c \ 



(21) 



1 Uy/C 

1 ujc 
\0 1 



(22) 



and 



Vy 

Vz 



k^ 



ky 
fcj. 



(23) 



Whence, the decomposition v = u + k can be viewed as a change in reference frames 
where an observer comoving with the volume element of the fluid whose hydrody- 
namic velocity is u will measure a given molecule's velocity as k while an observer 
in the laboratory sees the molecule moving at velocity v as given by Eq. (12 ip . 



III. RELATIVISTIC KINETIC THEORY 

In this section we will address the properties of a dilute, neutral, non-degenerate 
gas within the realms of special relativity. This system is thus described in a 
Minkowsky space-time whose metric is given by ds"^ = dx^ + dy'^ + dz^ — cdt^. For 
the molecules in this gas, the molecular four-velocity is given by 



V^ = 7^„ {W, C) 



where 



7u 



iM 



w 



2\ -1/2 



(24) 



(25) 



and w is the velocity. The distribution function has the same interpretation as 
above, being f {x'^ ,v^) d^xd^v the occupation number of a phase space cell. The 
special relativistic Boltzmann equation in the absence of external forces is given by 



vV,a = f = Aff) 

where the collision term is defined as 



J (/, f') = 1 1 if 'h'-fh) FadQdvl 



(26) 



(27) 



Here F is an invariant particle flux [7| which plays the role of the relative velocity, 
a is the collision cross section and the invariant differential volume in velocity space 
is dv* = ^. 

Here, as in our previous works, the proposed solution method for the kinetic 
equation is the Chapman- Enskog procedure to first order in the gradients [12]. As 



has been shown elsewhere |13J, this solution leads to a constitutive equation for the 
heat fiux in terms of gradients of the state variables. This is consistent with Onsager's 
regression of fiuctuations hypothesis and thus predicts no pathological behaviors in 



the system of hydrodynamic equations 



14| . Additionally, this system of equations 



to first order in the gradients, as predicted by kinetic theory, has been shown to 
present no causality issues both in the non-relativistic and relativistic cases [U]. In 
this case, the local equilibrium function is the Jiittner function 

f''' = ^VTY-pf^l (28) 

where lA^ = 7^ (m, c) is the hydrodynamic four- velocity, z = ksT /mc^ is the relativis- 
tic parameter and /C„ is the n-th order modified Bessel function of the second kind. 
As in the non-relativistic case, the transport equations are obtained by multiplying 
Eq. ( l26l) by collision invariants, in this case ip = 1, f^'. Indeed, the corresponding 
transport equation in the absence of external forces is 



v'^^fdv* 



(29) 



which yields the continuity equation for ^ = 1, the energy- momentum balance equa- 
tion for ip = mv^; that is, the momentum balance in the absence of external forces 
for /3 = 1, 2, 3 and the energy balance for /3 = 4. Equation ( l29l) can be expressed in a 
more conventional form as a general conservation law for four-fiows by defining the 
particle and stress-energy fiuxes as 

N" = [v'^fdv* (30) 



Tf"" = m v^v^fdv* (31) 

respectively. Thus, the transport equations are given by A^.^ = and T.^'' = 0. It is 
then appeahng to write the integrals in Eqs. fl5U]) and (IHT]) in terms of systematic 
and chaotic quantities in order to separate the different contribution to the fluxes, 
as done in the previous section (see Eqs. (IT5|) to ( !20|) ). To accomplish this, an 
appropriate transformation law has to be assigned in order to introduce the chaotic 
velocity. It is important to recall at this point that the hydrodynamic velocity is a 
local equilibrium quantity and is thus only defined in each differential volume where 
local equilibrium is assumed. If we fix our attention in a single random molecule, 
we can consider two reference frames, one in the laboratory (S) and one fixed in 
the volume where the molecule is contained. This second frame (S), in which the 
molecules would be seen static on the average, is moving with a speed u as seen by an 
observer fixed in S. Thus, observers in S and S would report that the corresponding 
velocities are given by 

f" = 7fc(^,c) (32) 

and 

orentz boost with velocity u and K"" = 7^ ( fc, c ) is 
We wish to remind the reader at this point that 
the contravariant transformation given in Eq. 033p is equivalent to the relativistic 
velocity addition law. 

With the transformation given in Eq. (1331) . Eqs. (130!) and (13T]) can be written as 



respectively. Here C^ is a 
the chaotic four- velocity |ll|. 



m = C^^ K'^fdK* (34) 

T^^ = m£^£^ I K'^K^fdK* (35) 



10 



where use has been made of the fact that, since dv* is an invariant quantity, dv* = 
dK*. Also the equihbrium distribution function given in Eq. fl28|) can be written in 
terms of the chaotic speed by use of the invariant 7^ = U^Vp/(? in a similar fashion as 
in the the non relativistic case where the argument of the Maxwellian is proportional 
to k'^. These two properties which allow the calculation of integrals in terms of K 
are verified in the Appendix. In order to obtain a general expression for T^'^ , we 
introduce an irreducible decomposition relative to the hydrodynamic four-velocity 
direction. That is, in this 3+1 representation a second rank tensor can be expressed 

as y 

rpt^y = rWW + T^'W + T^'W + r^'^ (36) 

where r^lA^j, = and r^^'^Uv = 0. The scalar, first and second rank tensors introduced 
can be expressed in terms of T'^" as 



r = T'^-^ (37) 

r^ = -^KT'^^Up (38) 



c^ 



r 



p,u 



h^h'^oT'^^ (39) 



respectively. Here h^'^ = g^^ + U^U^ j(? is the well known projector which satisfies 
U^Ji^ = 0. It is important to point out in this stage that in the phenomenological 
treatment, the quantities above are identified as the internal energy, heat flux and 
stress tensor without a kinetic theory based justification |5|. These definitions are 
in turn used in most kinetic treatments 2|, y, |6|. It is precisely the aim of this work 
to deduce, from purely kinetic grounds, that these quantities are indeed related to 
internal energy, heat flux and stress interpreted as averages over chaotic velocities in 
a similar fashion as in Eqs. (IT8l) - (l2(Ill . 
The scalar r can be calculated as 

c4 



m^^ / v^v'^fdv* = m I -flfdK* (40) 



11 



which is the internal energy per particle. To see that this is so, consider Eq. ( !29|) 
with ijj = mv'^ 

|- (m fv\^fdv*\ +-S^(m fvVfdvA = (41) 

where here, as in the rest of this work, latin indices run from 1 to 3 only. It is clear 
from Eq. (HT]) that the integral in the first term is indeed the total energy while the 
second integral is the energy flux. Thus, the equivalent to the total energy moment 
calculated in a rest frame yields the internal energy only, that is 

n€ = mc^ f -flfdK* (42) 

and thus, 

..-.„„(3..||) (43, 

For the vector quantity t'^ we have 

t" = -\k^p fv^yf^fdv* (44) 

which, using again the fact that W/jf '^ = —c^'jk can be expressed as an integral over 
the chaotic velocities as follows 

T'^ = h^^Cljj,K^fdK* (45) 

It can be shown (see the Appendix) that the contraction of the projector with the 
Lorentz transformation yields a tensor TZ^ = h'^C" given by 

7^^ = (46) 

K = C^, fora = l, 2, 3 (47) 



and thus 



^^j^kK^fdK* (48) 
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We now introduce the Chapman- Enskog expansion 

T^ = n^J ^kK^f^^'^dK* +%j jkK^f^'UK* (49) 

and notice that the first terms vanishes since, for /3 = 1,2,3 the integral 
TZ'i f 'jkK'^ f^^'dK* is odd in k and the /3 = 4 term in the sum is zero because 
7^4 = for any /i. Thus, only the integral with f^^^ survives and we can write 

T^' = n'; l^kK^f^'^dK* (50) 

In order to re-introduce the Lorentz transformation, we notice that 

j ^kK^f^^dK* = Q (51) 

since the internal energy, as all state variables, is obtained only through the equi- 
librium solution. That is, the subsiadiary condition, which the solution f^^^ will be 
enforced to satisfy, requires 

f -flf^^^dK* = for z 7^0 (52) 

Using Eq. (15T|) , one can write Eq. (|50|) back in terms of Ct which yields 

T'^ = C';J^kK''f^'UK* (53) 

By inspection of Eq. (HT]) one concludes that the integral J 'jkK^'f'^^^dK* is the heat 
flux in a rest frame where t>" = fC", and thus 

q^,^=c' j^uK^f^'^dK* (54) 

This expression is analogous to the one found in the non-relativistic case and full of 
physical content. The heat flux is physically the average flux of the chaotic energy, 
and Eq. (15^ is completely consistent which this idea. Now, in an arbitrary frame 

^' = ^2^Xo] (55) 

13 



which is, to the authors' knowledge, the first time that the heat flux is obtained only 
from a kinetic theory standpoint as the average of the peculiar kinetic energy flux of 
the molecules. 

For the second rank tensor in the stress-energy tensor decomposition, we calculate 
fromEq. 1^ 

rt"" = mWJffiC'^Cl f K^K^fdK* (56) 

or 

t"" = mTZ^n's I K'^K^ (/^^^ + f^^) dK* (57) 

For the local-equilibrium term we have 

m7^;;7^^ I K'^K^f^^dK* = mUX f K'^K'^f^^^dK* (58) 

Since f^^^ is even in k, only the a = b terms survive and thus 

m7^^7^^ f K'^K^f^^dK* = ph"'' (59) 

where we have introduced the well known result for the hydrostatic pressure 

p = m f {K')'f'UK* = m f [K^f f^'^^dK* = m f {K^f f^'^^dK* (60) 

and 

p = nksT (61) 

together with the identity 

n'^ni5''^ = c^cis""^ = h"" (62) 

Equation 
Weinberg 



was obtained in a similar fashion (using Lorentz transformations) by 
15] . nevertheless he did not address the dissipative case following a kinetic 



theory approach. 
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For the dissipative term in Eq. (157|) . which we write as U'^", we have 

W = mn'^Ws I K'^K^f^^dK* = mC^C"^ f K^K^f^^dK* (63) 

If Y&, is the Navier-Newton tensor calculated in a frame where the fluid is at rest 
nf,^ = m/i^/i^ f K^K^fdK* = mb^X f K^K^fdK* 

since in such frame h'^ = and h'^ = S^. Thus, the second rank tensor introduced in 
the stress-energy tensor is 

where 

n^'^ = r^r^njf (64) 

is the Navier-Newton tensor in an arbitrary frame. 

IV. SUMMARY AND FINAL REMARKS 

In the previous section, the different contributions to the stress-energy tensor 
for a single component, dilute gas in the framework of special relativity have been 
calculated by separating hydrodynamic and chaotic contributions to the molecular 
velocities. This has been accomplished by introducing Lorentz transformations to 
relate the velocity of a molecule as measured by an arbitrary observer with the one 
measured within a differential volume moving at the corresponding hydrodynamic 



velocity, an idea introduced by two of us in Ref. HJ, combined with Eckart's de- 
composition (5|. 

The main results of this work can be summarized in the fact that all quantities 
appearing in Eq. (136|) have been obtained strictly from kinetic theory using the 
concept of chaotic velocity. The first two terms are the equilibrium parts of the stress- 
energy tensor and are well known. The main accomplishment of the calculation here 

15 



shown are the dissipative terms which appear here in a natural way as averages over 
kinetic energy and momentum fluxes once the transformation between molecular and 
peculiar velocities is introduced. Also it has been shown that the heat flux transforms 
as a first rank tensor. 

A kinetic derivation of the stress-energy tensor for a dissipative fluid from first 
principles in kinetic theory has been lacking for some time and thus hindering a clear 
derivation of the relativistic Navier-Stokes equations. Equations ( 1M|) and (|35|) satisfy 
both needs and, in turn, pose a new question. Since both heat and momentum fluxes 
in Eq. ( l36l) are given by Eqs. ( l55l) and ( 1M|) respectively, the hydrodynamic velocity 
factors introduced by the Lorentz transformations will induce new non-linearities in 
the system of hydrodynamic equations. This could yield new relativistic effects for 
the relativistic gas which may be measurable. This question and will be addressed 
in the future. 

Appendix 

In this appendix the relations 



dv* = dK* (65) 



and 



H'v, = 7, (66) 

are shown to hold where W^, v^ and K^^ are the hydrodynamic, molecular and chaotic 
four-velocities respectively. Also, we verify that the tensor quantity TZ'^^ is indeed 
given by Eqs. (jl6]) and (H7|) . 

To verify Eqs. (165|) and (166|) . we consider two reference frames S and S with 
a relative speed u with respect to each other. That is, S may be considered the 
laboratory frame while 5 is a frame fixed to a volume element in the fluid. For the 

16 



sake of simplicity, we take the x direction parallel to u. In this situation, we have 
three four- vectors related to a given molecule in such fluid element 

K" = 7fc ( A;, c ) velocity of the molecule as measured by an observer in S 

v'^ = 7t, {w, c) velocity of the molecule as measured by an observer in S 

U^ = 7„ {u, 0, 0, c) relative velocity between S and S 

The relationship between tensors in both references frames given by the Lorentz 
transformation 



(67) 



then 





7« 


c '« 


M _ 


1 







f7« 1 


lu 



A^ = C^A" (68) 

Since the molecule's velocity, as measured in S, is v'^ = K" we have 



v^ = C^K" = 7 J 7„ (m + ki) , k2, /cs, 7«c 1 + 



uki 



(69) 



In order to show the invariance of the volume element dv* = cd^v/v^ we start 
from 



where the Jacobian is given by 



d'^v = Jd-^K 



J = det 



Qya 



(70) 



(71) 
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and is calculated as follows 

Qya Q 



d 



QKb ORb f^"^'] ^'"dK" f^' 




a = 1 
a^ 1 



where use has been made of the fact that, since K'^ is a four- velocity, K^^K^ = —(? 
and thus 






Q^b ^QKb "OK^ ^dK^ 



Then, the Jacobian is 



J = 7.(l-^§)=^7«(i^4-^i^i)=^7.7.(c+^/ci)=.^ 



and thus 






Regarding the scalar product W^Vi, we have 

which can be readily calculated using that 

Uv = lu {u, 0, 0, -c) 
and the transformation given in Eq. (167|) as follows 



U^v'^ = juuCiK^ - c^uCtK^ = ll 



u 



c 



7' 



U 



-1 



and thus 



U^v" = -cK' 



(72) 

(73) 
(74) 

(75) 
(76) 

(77) 
(78) 



The results in Eqs. (17^ and (178|) allow the calculation of moments of the distribution 
function in terms of the chaotic velocity in as similar way as in the non-relativistic 
case: 

(79) 



exp(^)^rf.* = |exp(-|):7rfir* 
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where J is an arbitrary tensor. 

Now we turn to the proof of Eqs. f H6|) and (H7|) . Firstly, since £4 = ~) 

K = K'^t = /i^— = (80) 

To obtain Eq. (HTj), we separate two cases. For /i = 4, since h\ = l — 'y'^ and hf = 7-^ 

K = htc: = {i-^^)^+^^ct (81) 

For the second term we use that 

^b ^b ^b I rb , ^a^ \ ^a 



and thus 



^£^ = ^ U^ + ^7^^ = 7- (82) 

c c V c2(7 + 1)7 c ^ ^ 



^a = ^ = ^^ (83) 



Finally, for /x = £ = 1, 2, 3 

7^^ = hie: = hiCl + /^i£^ (84) 

or, using that h^ = —7— 

^a = y^l + ^^ J ^a - 1-^ (85) 

Now, by introducing Eq. 082|) in Eq. (ISSj) one obtains 

7^1 = ^'a (86) 

This completes the proof. 



[1] F. Jiittner; Ann. Physik und Chemie 34 (1911) 856. 
[2] W. Israel; J. Math. Phys 4, 1163 (1963). 

19 



[3] S. R. de Groot, P. Mazur; Non-Equilibrium Thermodynamics, Dover Publications NY 
(1984). 



[4 
[5 
[6] 

[7] 



[9 
[10 

[11 
[12 

[13 

[14 

[15 



S. Chapman, T. G. Cowling; The Mathematical Theory of Non-Uniform Gases. 

C. Eckart, Phys. Rev. 58, 267 (1940); ibid 58, 919 (1940). 

S. R. de Groot, W. A. van Leeuwen and Ch. van der Wert; Relativistic Kinetic Theort, 

North Holland Publ. Co., Amsterdam (1980). 

C. Cercignani and G. Medeiros Kremer; The Relativistic Boltzmann Equation: Theory 

and Applications, Cambridge University Press 3rd Ed., UK (1991). 

R. Calusius, Annalen der Physik 100, 353-380 (1857). 

S. Brush; The kind of Motion we call Heat, North Holland, Amsterdam (1986) . 

J. C. Maxwell, On the dynamical theory of gases, in Scientific Papers of J. C. Maxell; 

W. D. Niven, editor Vol. 2 p. 26 Dover Publications NY (1965). 

A. Sandoval- Villalbazo, L. S. Garcia-Colin, Physica A 278 428 (2000). 

A.L Garcia-Perciante, L. S. Garcia-Colin and A. Sandoval- Villalbazo, Gen. Rel. Grav. 

41 1645 (2009). 

A. Sandoval- Villalbazo, A. L. Garcia-Perciante, L. S. Garcia-Colin, Physica A 388, 

3765 (2009). 

A. L. Garcia-Perciante, A. Sandoval- Villalbazo, Journ. Non-Newtoniana Fluid Mech. 

available online March 2010. 

S. Weinberg; Gravitation and Cosmology, Wiley & Sons, USA (1972); Chapter 2. 



20 



